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Abstract 

We study the coherent atomic tunneling between two zero-temperature Bose- 
Einstein condensates (BEC) confined in a double- well magnetic trap. Two 
Gross-Pitaevskii equations for the self- interacting BEC amplitudes, coupled 
by a transfer matrix element, describe the dynamics in terms of the inter-well 
phase-difference and population imbalance. In addition to the anharmonic 
generalization of the familiar ac Josephson effect and plasma oscillations oc- 
curring in superconductor junctions, the non-linear BEC tunneling dynamics 
sustains a self-maintained population imbalance: a novel "macroscopic quan- 
tum self-trapping effect" . 

PACS numbers: 03.75.Fi,74.50.+r,05.30.Jp, 32.80.Pj 
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The recent experimental observation of the Bose-Einstein condensation (BEC) in a dilute 
gas of trapped atoms , has generated much interest in the properties of this new state 
of matter. A quite fascinating possibility is the observation of new quantum phenomena 
on macroscopic scales, related with the superffuid nature of the condensate. In fact, bro- 
ken symmetry arguments show that the condensate atoms can be described by a common, 
"macroscopic" one-body wave function ty(f,t) = Jp~ e %e (the order parameter), with p the 
condensate density. For a weakly interacting BEC, the order parameter obeys a non-linear 
Schroedinger, or Gross-Pitaevskii equation (GPE) ||: 

where V ext is the external potential and g = 47T ^ a is the inter-atomic scattering pseudo 
potential, with a, m the atomic scattering length and mass respectively. 

The GPE has been successfully applied to investigate the collective mode frequencies 
of a trapped BEC in the linear regime the relaxation times of monopolar oscillations 
H, and, because of the non-linear self-interaction, it could also induce chaotic behavior in 
dynamical quantum observables [^,0 

The existence of a macroscopic quantum phase (difference) was dramatically demon- 
strated recently [[|. A far off- resonant intense laser sheet divided a trapped condensate, 
creating a large barrier in between. Switching off the double-well trap, the two released 
condensates overlapped, producing a robust "two-slit" atomic interference pattern, clear 
signature of phase coherence over a macroscopic scale (~ 10~ 2 cm). The non destructive 
detection of phase-differences between two trapped BEC could be achieved by lowering the 
intensity of the laser sheet. This allows for atomic tunneling through the barrier, and the 
detection of Josephson-like current-phase effects |pJ,0-^[. In superconductor Josephson junc- 
tions (SJJ), phase coherence signatures include a dc external voltage producing an ac current, 



or the "plasma" oscillations of an initial charge imbalance [10,111. For neutral superffuid 



He II, voltage drives, tunnel junctions, or capacitive charges are absent. The only accessi- 



ble Josephson analogue |12[ involves two He II baths connected by a sub-micron orifice, at 



which vortex phase-slips JT3 support a chemical potential (height) difference, through the 
Josephson frequency relation. 

Although the trapped BEC is also a neutral-atom Bose system, its population can be 
monitored by phase-contrast microscopy; the double-well curvatures and barrier heights can 
be tailored by the position and the intensity of the laser sheet partitioning the magnetic 
trap f2j. We note that the chemical potential between the two condensates depends both 
on the zero-point energy difference, that acts like an external "dc" SJJ voltage, and on the 
non-linear interaction that, through an initial population imbalance, acts like a capacitive 
SJJ charging energy. Thus, we suggest that the BEC tunnel junction can show the analogues 
of the familiar Josephson effects in supercunductor junctions, with the ability to tailor traps 
and the atomic interaction compensating for electrical neutrality. 

In this Letter we study the atomic tunneling at zero temperature between two 
non-ideal, weakly linked BEC in a (possibly asymmetric) double well trap. This in- 
duces a coherent, oscillating flux of atoms between wells, that is a signature of the 
macroscopic superposition of states in which the condensates evolves. The dynamics is 
governed by two Gross-Pitaevskii equations for the BEC amplitudes, coupled by a transfer 
matrix element (Josephson tunneling term). Analogues of the superconductor Josephson 
effects such as the ac effects and "plasma" oscillations are predicted. We also find that the 
non-linearity of the dynamic tunneling equations produces a novel self-trapping effect. 

Consider a double- well magnetic trap 1,2 as in Fig.(l). This system can be described by 
a two-state model: 

ih^ = {E\ + C^JVOVi - Kfo (2a) 
tk^ = (E° 2 + U 2 N 2 )^ 2 - (26) 



with uniform amplitudes ipi j2 = yA^i,2 e i6>1 ' 2 , where iV lj2 , 1)2 are the number of particles 
and phases in the trap 1,2 respectively, and K is the coupling matrix element jnj. The 
parameters E® 2 , ^1,2 and K can be determined from appropriate overlap integrals of the 
time-independent GPE eigenfunctions $12 of the isolated traps, as outlined later. The total 
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number of atoms Nt = N\ + N 2 is constant, but we stress that a coherent phase description, 
i.e. the existence of definite phases # 12 , implies that the phase fluctuations (~ i — WM) 
must be small, giving Aq j2 > N min ~ 10 3 , say. 

We note that the BEC inter-trap tunneling Eq.s(2) are similar in form to models of 
single-trap atomic- level transitions ||, or polaron hopping in semiclassical approximation 
16fl , although they describe quite different physics. 

In terms of the phase-difference <fi = 9 2 — Q\ and fractional population difference — 1 < 
z = < i ; Eq.s(2) become (h = 1): 



z = — Vl — z 2 siri(f) (3a) 



Az + cosch + AE (36) 

:i-z*) 



where the time has been rescaled as 2Kt — > t. The dimensionless parameters are: 

AE = {El - El)/{2K) + {U x - U 2 )N T /{AK) (46) 

A=(U 1 + U 2 )N T /{AK) (4c) 

For two symmetric traps, = E% (AE = 0), Ui = U 2 = U, and A = UN T /2K. In the 
following we will assume a positive scattering length a (A > 0); note, however, that Eq.s(3) 
are invariant under the trasformation A — > —A, — > —0 + 7r, AE AE. 

The z, <fi variables are canonically conjugate, with i = — = and the Hamiltonian 
is given by: 

A 



H = —z 2 — Vl-z 2 cost + AEz (5) 
In a simple mechanical analogy, H describes a non — rigid pendulum, of tilt angle <fi and 



length proportional to vl — z 2 , that decreases with the "angular momentum" z. 
The inter-trap tunneling current is given by: 

I = z-L = I^(l- z z)sin(i>- I = KN T (6) 

The detailed analysis of Eq.s(3) with exact analytical solutions will be presented elsewhere 



17]]; here we consider three regimes. 



1) Non — interacting limit. For symmetric wells and negligible interatomic interactions 
(A — > 0), Eq.s(2) can be solved exactly, yielding Rabi-like oscillations in the population of 
each trap with a frequency: 

Ur = 2K (7) 

as studied in |7]||. However, the ideal Bose gas limit is not accessible experimentally. 

2) Linear regime. In the linear limit (\z\ << 1, \<p\ << 1) Eq.s(3) become: 

i ~ -(j) (8) 

</>~(A + l)z (9) 

These describe the small amplitude oscillations of the pendulum analogue, with a sinusoidal 
z(t) with a frequency: 



u L = \J 2UN T K + 4K 2 (10) 

The BEC oscillations of population should show up as temporal oscillations of phase-contrast 
patterns 0, or other probes of atomic population [|TJ. 
Linearizing Eq.s(3) in z(t) only, we obtain: 

z ~ —sincj) (11°) 

(j) ~ AE + (A + cos<f))z (116) 
I ~ I sin(f) (He) 

For large trap asymmetries with AE » (A + cos(<ft))z, we have: <fi = 0(0) + AE't, giving 
an oscillating z(t) with frequency: 

u ac ~ ^ - £ 2 ° (12) 

where an "ac" current /(£) is produced by the "dc" trap asymmetry AE. It is simple 
to show that a small oscillation in the laser position, or in its intensity (K — > K(l + 
5sin(uot)), 5 << 1), will result in a dc inter-trap current of non-zero time average < 
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I(t) >~ 5 < sin(uJot)sin(uj ac t) >^ 0, at resonant match uoq = 0J ac . This is the analogue of 
the Shapiro effect |10[ in superconductor junctions. 

In SJJ, the current of Cooper pairs iVi )2 is Isjj = —2e(Ni — N 2 ) = 2eEj sincp, and 
the Josephson frequency relation for the relative phase is <fi = Afi = 2eV + (Aq — N2)E C , 
for a dc applied voltage V and a junction capacitance C, with E c = (2e) 2 /2C |TIJ. These 
SJJ equations can be directly compared with the BEC Eq.s(3). It is then clear that the ac 
Josephson frequency uj ac = 2eV and the Josephson "plasma" frequency uj p = y/E c Ej 



are the analogue of Eq.(12) and Eq.(10) respectively. Note however that, for SJJ, u p is 
independent of the barrier cross section A, since Ej ~ A and E c ~ C~ l ~ A^ 1 , while the 
BEC has ul — A 1 / 2 since K ~ A, and the bulk energy UNt is approximately A- independent. 

3)Non— linear regime. A numerical solution of Eq.s(3) yields non-sinusoidal oscillations, 
that are the anharmonic generalization of the Josephson effects. Moreover, an additional 
novel non-linear effect occurs in the BEC: a self-locked population imbalance. 

Fig. (2) shows solutions of Eq.s(3) with initial conditions z(0) = 0.6, (f)(0) = and 
illustrative parameters A = 1, 8, 9.99, 10, 11 respectively. The sinusoidal oscillations around 
z = became anharmonic as A increases, Fig.(2a,b), and with a precursor slowing down, 
Fig. (2c), there is a critical transition for A = A c = 10, dashed line in Fig(2d). Then for 
A = 11 the population in each trap oscillates around a non-zero time averaged < z(t) >^ 0, 
solid line in Fig. (2d). In the non-rigid pendulum analogy, this corresponds to an initial 
angular momentum z(0) sufficiently large to swing the pendulum bob over the <fi = ir vertical 
orientation, with a non-zero < z(t) > average angular momentum corresponding to the 
rotatory motion. This critical behavior depends on A c = A c [2;(0), 0(0)], as can be easily 
found from the energy conservation constraint. In fact from Eq.(5), the value z(t) = is 
inaccessible at any time if: 

Jl - z(0) 2 cos(6(0)) + 1 
A > 2 (V ) (13 ) 

The full dynamical behavior of Eq.(3) is summarized in Fig(3), that shows the z — <fi phase 
portrait with constant energy lines. There are energy minima along z = at 2n7c, and 
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"running" solutions < <p(t) >y^ with < z(t) >y^ 0, moving along the sides of these wells. 
The vertical points (p = (2n + l)ir, that would be isolated unstable points for the rigid 
pendulum, now support oscillations of restricted range, as a consequence of non-rigidity, i.e. 
non-linearity. 

The self-trapping of an initial BEC population imbalance, seen in Fig. (2d) and Fig. (3) 
arises because of the inter-atomic interaction in the Bose gas (non-linear self-interaction 
in GPE). It has a quantum nature, involving the coherence of a macroscopic number of 



atoms. It differs from single polaron trapping of an electron in a medium and from 



gravitational effects on He II baths [fHql . It can be considered as a novel "macroscopic 
quantum self-trapping" (MQST). 

Non-linear effects like MQST are unobservable in SJJ. The requirement that the chemical 
potential difference fix — fi 2 — (Ni — N 2 )E C must be less than the quasiparticle gap 2A qp (to 
prevent a jump off the Josephson I-V branch) implies a stringent constraint in imbalances 
\z\ < (2 A qp / E c Nt) — 1CT 12 for typical parameters. For the BEC, the requirement that 
tunneling does not access excitation energies is much less restrictive. As an example, let 
consider two weakly linked condensates of Nt — 10 4 atoms, confined in two symmetric 
spherical traps with frequency ujq ~ 100 Hz. Evaluating Eq.s(17) below with a simple 
variational wave function jrj|, we have E° ~ 0.5 nK, UN? — 3 nK, and from an estimation 
of the excitation gap we obtain the constraint \z\ < 0.5. Taking K ~ 0.1 nK, we have 
A ~ 10, close to the onset of the critical behavior. Increasing the number of particles 
or the trap frequencies, the value of A can increase, making MQST observable. Typical 
frequency oscillations are ujl — KHz for Nt — 10 6 , that should be compared with the 
plasma frequency of SJJ that are 0] of the order of u p ~ GHz. 

We now outline the derivation of the parameters Eq, U, K. 

To this purpose we note that in the barrier region the modulus of the order parameter 
in the GPE is exponentially small. This allows us to look for an eigenstate of GPE of the 
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form: 

* = Mt)*i(x)+Mt)*2(?) (14) 

where $1, $2 are the ground state solutions for isolated traps with N\ = iV 2 = Nt/2. 
Replacing Eq.(14) in the GPE Eq.(l) with the conditions: 

J $i$ 2 df ~ (15) 

J \\$ 1 \\ 2 dr = J ||$ 2 || 2 dr=l (16) 
r h 2 

£° 2 = J [^IV^I 2 + ^l 2 Ve X t]df (17a) 

C/1,2 = So / $i, 2 df (176) 

^ = - / [^(V$iV$ 2 ) + $^2^]^ (17c) 



and 



we obtain: 



At finite temperature the interaction between the normal component of the Bose gas 
with the condensate should be included, and the parameters become temperature dependent. 
Such corrections are small for temperatures smaller than excitation energies. High density 
BEC could induce quasiparticle/collective mode scattering with finite lifetime of the coherent 
oscillations ||; phase diffusion could induce phase coherence collapse and revival fL9| . These 
effects deserve further studies. 

In conclusion, the BEC coherent atomic tunneling in a double-well trap induces non- 
linear population oscillations that are a generalization of the sinusoidal Josephson effects 
familiar in superconductors. A novel population imbalance occurs for parameters beyond 
critical values: a macroscopic quantum self-trapping effect. 

Discussions with S.Raghavan and useful references from L. Bonci and G. Williams are 
acknowledged. 
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FIGURES 

FIG. 1. The double well trap for two Bose-Einstein condensates with Nip and E®p the number 
of particles and the zero-point energies in the trap 1,2 respectively. 

FIG. 2. Fractional population imbalance z(t) versus rescaled time, with initial conditions 
z(0) = 0.6, phase difference 0(0) = 0, and A = 1 (a), A = 8 (b), A = 9.99 (c), A = 10 (dashed line, 
d), A = 11 (solid line, d). 

FIG. 3. Constant energy lines in a phase-space plot of population imbalance z versus 
phase difference 0. Bold solid-line: z(0) = 0.6, 0(0) = 0, A = 1,8,10,11,20. Solid line: 
z(0) = 0.6, 0(0) = 7T, A = 0, 1, 1.2, 1.5, 2. 
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